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Abstract.  Turbulent  channel  flow  with  a rib- roughened  wall  has  been  com- 
puted with  constant  heating  rate  using  large  eddy  simulation.  An  implicit, 
second  order  accurate  finite  volume  scheme  was  used  to  solve  the  time 
dependent  filtered  set  of  equations  to  determine  large  eddies,  while  a dy- 
namic subgrid-scale  model  was  used  to  account  for  the  subgrid  scale  effects. 
A dynamic  Prandtl  number  was  introduced  to  calculate  subgrid  scale  heat 
transfer.  The  effects  of  strong  heating  on  the  rib-roughened  wall  were  in- 
vestigated and  Nusselt  numbers  were  computed  and  compared  with  exper- 
imental data.  The  instantaneous  plots  revealed  the  influence  of  turbulent 
heat  transfer  on  velocity,  temperature  and  vorticity  in  this  particular  ge- 
ometry. 


1.  Introduction 

Engineers  are  using  three  different  approaches  to  cope  with  turbulence.  In 
industry,  Reynold  average  Navier-Stokes  (RANS)  methods  employing  two 
equation  models  are  widely  used  in  a variety  of  forms.  Direct  numerical 
simulation  is  mainly  used  in  academic  research.  Large  Eddy  Simulation 
(LES)  which  originated  in  academia  tends  to  be  a more  and  more  useful 
analysis  tool  for  industry  in  recent  years.  Several  commercial  codes  (such  as 
FLUENT  and  StarCD)  have  implemented  some  types  of  LES  models.  LES 
resolves  all  the  energy-carrying  large  eddies,  but  models  small  eddies  with 
length  scales  smaller  than  the  grid  size.  This  equips  LES  with  the  ability 
to  analyze  complex  three  dimensional  time  dependent  flows  at  a moderate 
computational  cost  where  RANS  simulations  usually  fail. 
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The  rib-roughened  channel  flow  is  a good  test  case  for  LES  simulation 
due  to  its  complex  flow  features,  such  as  flow  separations,  recirculation 
zones,  and  re-attachment.  Cooling  in  gas  turbine  blades  is  a potential  ap- 
plication for  this  kind  of  flow.  Early  work  was  done  by  Ciofalo  and  Collins 
[1]  for  a rib- roughened  channel  using  the  Smagorisky  model  and  a commer- 
cial CFD  code  employing  an  incompressible  flow  formulation.  Yang  and 
Ferziger  [2]  used  a dynamic  subgrid-scale  model  to  simulate  a turbulent 
flow  over  an  obstacle  the  geometry  of  which  is  different  from  Ciofalo’s  case. 
Dailey  [3]  developed  a compressible  formulation  for  large  eddy  simulation 
and  applied  it  to  the  simulation  of  a compressible  channel  flow  with  heat 
transfer.  Further,  Meng  and  Pletcher  [4]  calculated  an  isothermal  flow  case 
similar  to  that  of  Ciofalo  and  Collins  but  using  dynamic  SGS  model  and 
a finer  grid.  This  paper  is  a continuation  of  that  research  and  considers 
effects  of  heat  transfer  and  property  variations. 


2.  Problem  Description 


The  problem  of  interest  is  the  turbulent  flow  and  heat  transfer  in  a two- 
dimensional  plane  channel  with  a periodic  array  of  transverse  square  ribs 
on  one  wall.  The  simulation  was  set  up  such  that  comparisons  could  be 
made  to  the  experimental  data  of  Drain  et  al.  [5]  and  Bates  et  al.  [6],  who 
made  detailed  mean  and  fluctuating  velocity  measurements.  The  blockage 
ratio,  which  is  the  ratio  of  the  rib  height,  h,  to  channel  height,  H,  was 
h/H  = 0.2.  The  pitch-to-height  ratio  of  the  ribs  was  7.2,  where  the  pitch, 

P,  is  the  streamwise  spacing  of  the  ribs.  This  value  of  pitch-to-height  ratio 
was  consistently  found  empirically  to  yield  the  greatest  enhancement  of 
heat  transfer  rates  [1].  The  nominal  Reynolds  number  based  on  hydraulic 
diameter  (Dh  = AS)  and  bulk  velocity  was  Rep  = 20,000. 

For  simulations  with  heat  transfer,  a constant  heat  flux,  qw  = ^ ™ — , 

’ '-'pPrefUreflref  ’ 

was  applied  to  both  the  lower  wall  and  the  upper  wall.  On  the  front,  top, 
and  back  surfaces  of  the  rib,  a constant  heat  flux,  qw,rib,  was  applied  with 
a magnitude  that  was  one-third  of  the  heat  flux  on  the  lower  wall,  qw.  Con- 
sequently, the  total  heat  added  to  the  ribbed  channel  was  the  same  as  for 
a smooth  channel  with  no  ribs  and  wall  heat  flux,  qw.  For  the  simulations 
presented  here,  qw  = 0.002,  and  qw<rib  = 0.000667. 

The  computational  domain  only  contained  one  rib.  This  section  was 
assumed  to  be  embedded  in  a periodic  array  of  many  ribs  in  a fully  devel- 
oped flow,  which  permitted  the  use  of  periodic  or  step  periodic  boundary 
conditions  in  the  streamwise  direction. 
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2.1.  NUMERICAL  PROCEDURE 


A system  of  Favre-filtered  Navier-Stokes  equations  was  solved  using  a dy- 
namic sub-grid  scale  model.  The  details  can  be  found  in  Dailey’s  Ph.D. 
thesis  [3].  The  simulations  were  primarily  run  with  the  dual  time-stepping 
LU-SGS  scheme  with  the  solution  domain  divided  into  8 blocks  for  the 
coarse  grid,  and  16  blocks  for  the  fine  grid.  The  fine  grid  simulations  were 
run  with  17  processors  (one  for  master  node,  sixteen  for  slave  nodes)  on  the 
IBM-SP2  and  required  about  10  - 12  hours  of  wall  clock  time  per  10,000 
time  steps,  depending  on  the  number  of  subiterations  required  for  each  step 
in  physical  time. 

For  the  solid  walls,  including  the  upper  and  lower  walls  and  rib  surfaces, 
the  no-slip  velocity  and  zero  normal  pressure  gradient  boundary  conditions 
were  enforced.  The  wall  temperature  was  determined  from  the  specified 
heat  flux.  All  solution  variables  were  assumed  to  be  periodic  in  the  z di- 
rection, the  only  homogeneous  direction  for  this  flow.  In  the  streamwise 
direction,  the  ^-momentum,  v and  w velocities,  and  periodic  component  of 
the  pressure,  pp,  were  assumed  to  be  periodic.  The  mean  pressure  gradient 
parameter  was  adjusted  at  each  time  step  to  maintain  the  desired  mass 
flow  rate.  The  temperature  was  assumed  to  be  stepwise  periodic,  with  A Tx 
given  by 


A Tx 


QwLx 

2 (m/LyLz) 


(1) 


where  Ly  = 2 and  Lz  = n. 

After  the  turbulent  flow  was  deemed  be  statistically  stationary,  the  sim- 
ulation was  run  another  10,000  time  steps  to  compute  the  ensemble  aver- 
aged turbulence  statistics. 


2.2.  GRID  AND  RUN  PAR  AMETERS 

The  dimensions  of  the  computational  domain  were  7.2 h x 5 h x 7.85 h,  or 
2.885  x 2<5  x 7t<5,  in  the  x.  y,  z directions,  respectively.  The  coarse  grid  used 
40  x 32  x 24  control  volumes,  with  8 cells  on  each  surface  of  the  rib.  The 
fine  grid  consisted  of  80  x 64  x 48  control  volumes,  with  16  cells  on  each 
surface  of  the  rib.  The  fine  grid  is  shown  in  Fig.l.  The  nondimensional  time 
step  was  0.01. 


3.  Results 

3.1.  MEAN  AND  FLUCTUATING  VELOCITY  PROFILES 

The  mean  streamwise  velocity  profiles,  normalized  by  the  average  bulk  ve- 
locity on  the  top  of  the  rib,  are  depicted  in  Figs.  2 and  3 at  planes  A and 
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B,  respectively.  Plane  A was  located  mid-way  between  the  ribs,  and  plane 
B was  located  on  the  middle  of  the  rib  top  surface.  Both  experimental  data 
and  current  simulation  results  are  shown  to  facilitate  the  comparison.  From 
these  two  figures,  it  can  be  seen  that  both  coarse-grid  results  and  fine-grid 
results  agree  reasonably  well  with  the  experimental  data  in  predicting  av- 
eraged velocity,  although  fine-grid  results  are  slightly  better. 

The  root-mean-square  ( rms ) of  the  streamwise  velocity  fluctuations  are 
shown  in  Figs.  4 and  5.  In  the  prediction  of  velocity  fluctuations  in  section 
A,  the  fine-grid  result  shows  much  better  agreement  with  experimental 
data,  even  though  both  results  underestimate  the  velocity  fluctuations  to 
some  extent.  For  section  B,  the  fine-grid  result  captures  precisely  the  peak 
value  near  the  top  of  the  rib,  and  continues  to  be  closer  to  the  experimental 
data  than  the  coarse-grid  result.  It  seems  that  the  fine  grid  is  necessary  to 
calculate  the  turbulent  statistics,  e.g.  velocity  fluctuations,  accurately. 

3.2.  HEAT  TRANSFER  PARAMETERS 

The  local  Nusselt  number  distribution  is  shown  in  Fig.  6,  compared  to 
the  experimental  data  of  Lockett  and  Collins  [7]  at  two  different  Reynolds 
numbers.  The  distance  along  the  lower  wall  of  the  channel,  including  the 
rib  surfaces,  is  represented  by  x,  and  the  location  of  the  lower-front  corner 
of  the  rib  is  xle.  The  normalized  Nusselt  number  distribution,  as  shown  in 
Fig.  7,  used  the  average  Nusselt  number  as  the  normalization  factor,  which 
is  defined  as  below 


1 fs 

NuDiavg  = — / Nuo{s)ds 
& Jo 


(2) 


where  S is  the  distance  along  the  lower  wall,  including  the  rib  surfaces.  The 
LES  prediction  of  Nusselt  number  agrees  quite  well  with  experimental  data. 
It  reproduced  the  overall  streamwise  profile  of  Nup,  with  a local  maximum 
immediately  upstream  of  the  rib,  a global  maximum  immediately  adjacent 
to  the  left-upper  corner  of  the  rib,  and  a local  minimum  immediately  down- 
stream of  the  rib.  The  fine  grid  results  performed  slightly  better  than  coarse 
grid  ones,  in  particular,  in  the  prediction  of  the  peak  values.  In  Fig.  7,  use 
of  the  normalized  Nusselt  number  permits  comparison  of  results  obtained 
from  different  Reynolds  numbers.  The  examination  of  the  profiles  reveals 
that  the  LES  simulation  is  in  very  good  agreement  with  experimental  data, 
except  that  the  maximum  value  of  Nup  is  slightly  underestimated.  Note 
that  there  is  an  increase  of  Nusselt  number  near  the  trailing  edge  of  the  rib 
for  the  experimental  results  {Rep  — 29,870),  probably  due  to  a recircula- 
tion zone  that  develops  on  top  of  the  rib  at  high  Reynolds  numbers. 
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3.3.  ENSEMBLE  AVERAGED  AND  INSTANTANEOUS  CONTOURS 

Some  ensemble  averaged  and  instantaneous  contours  for  a rib-roughened 
channel  flow  with  heat  transfer  arc  presented  in  this  section. 

Figure  8 shows  mean  velocity  vectors  with  a ensemble  averaged  con- 
tour of  u velocity  as  the  background.  This  plot  was  created  based  on  the 
simulation  results  without  heat  transfer  for  the  purpose  of  comparison.  A 
large  recirculation  zone  behind  the  rib  and  a smaller  recirculation  zone  in 
front  of  the  rib  are  evident.  A large  fraction  of  the  region  between  the  ribs 
is  in  the  recirculation  zone  for  the  current  pitch-to-height  ratio.  Drain  and 
Martin  [5]  reported  a reattachment  length  of  4.3h.  The  present  simulation 
gives  a reattachment  length  of  4.2  — 4.3/i. 

A 3D  geometric  representation  of  this  flow  (instantaneous  u velocity) 
is  shown  in  Fig.  9.  A cutting  plane  is  located  midway  in  the  spanwise 
direction.  The  instantaneous  velocity  contours  are  very  irregular  because 
of  the  nature  of  turbulence. 

The  following  figures  show  the  contours  of  a variety  of  variables  on  two 
cutting  planes,  which  are  located  at  8%  and  92%  channel  height,  respec- 
tively, and  parallel  to  the  top  or  bottom  walls.  Figure  10  reports  the  in- 
stantaneous u velocity.  On  the  top  cutting  plane,  there  are  streak  patterns, 
which  are  characteristics  for  the  small-scale  turbulent  structure  near  the 
wall.  On  the  bottom  cutting  plane,  no  obvious  streak  pattern  can  be  found 
even  though  there  are  still  some  variations  of  u velocity.  The  existence  of 
recirculation  area  before  and  after  the  rib  is  thought  to  be  the  reason  for 
this  streak  deterioration. 

The  instantaneous  temperature  is  presented  in  Fig.  11.  The  top  cutting 
plane  is  uniformly  and  moderately  heated.  The  bottom  cutting  plane,  how- 
ever, has  a large  cool  region  upstream  of  the  rib  and  a small  hot  region 
downstream  of  the  rib  followed  by  another  cool  region. 

The  magnitude  of  vorticity  is  reported  in  Fig.  12.  There  exist  some 
streak  structures  of  vorticity  on  the  top  cutting  plane,  but  the  vorticity 
is  weak  and  smooth  on  the  bottom  cutting  plane  except  a narrow  area 
immediately  before  the  front  surface  of  the  rib.  The  recirculation  of  flow  is 
not  conducive  to  the  formation  of  streak  structure  of  vorticity. 

4.  Conclusions 

The  simulation  of  turbulent  flow  and  heat  transfer  of  air  in  a planar  channel 
with  transverse  square  ribs  on  one  wall  was  attempted.  Very  good  agree- 
ment between  fine  grid  results  and  experimental  data  were  obtained  in 
predicting  velocity  and  velocity  fluctuations.  Coarse  grid  results  agree  well 
with  experimental  data  in  velocity  profiles  but  not  as  well  as  the  fine  grid 
results  in  velocity  fluctuation  profiles.  The  Nusselt  number  profiles  obtained 
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from  the  current  simulation  agree  reasonably  well  with  experimental  data 
which  tends  to  validate  the  current  SGS  model  for  turbulent  heat  transfer. 
The  instantaneous  contours  of  some  important  variables,  e.g.  u velocity, 
temperature,  and  vorticity,  were  presented  and  analyzed.  It  is  found  that 
the  recirculation  of  flow  has  a negative  impact  on  the  formation  of  streak 
structure  near  the  wall. 
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Figure  1.  Computational  grid  for 
rib-rougliened  channel  flow 


Figure  2.  Mean  streamwise  velocity  nor- 
malized by  bulk  velocity  for  Section  A 


Figure  4-  rms  streamwise  velocity  fluctu- 
ations normalized  by  bulk  velocity  for  Sec- 
tion A 


Figure  5.  rms  streamwise  velocity  fluctu- 
ations normalized  by  bulk  velocity  for  Sec- 
tion B 
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Figure  3.  Mean  streamwise  velocity  nor- 
malized by  bulk  velocity  for  Section  B 


Figure  6.  Local  Nusselt  number  compared 
to  experimental  data  of  Lockett  and  Collins 
(1990) 
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Figure  7.  Normalized  Nusselt  number 
compared  to  experimental  data  of  Lockett 
and  Collins  (1990) 


Figure  10.  Instantaneous  u velocity 


